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Abstract 

Four sets of data determining parameters of ao(1450) and ao(980) are re-analysed. These are Crystal Barrel 
data for pp annihilation at rest to r]TY n in (i) liquid hydrogen and (ii) gas, and to (iii) K^-k^K^ and (iv) 
u-r + -r~-k° (mostly upn). Dispersive corrections due to opening of inelastic thresholds are treated fully. This 
stabilises parameters of ao(1450) substantially. The mass of its peak is 1448± 13(stat) ±25(syst) MeV and its 
mean full width at half maximum is 192±9±9 MeV. The pole position is M-iT/2 = 1432±13±25-i(98±5±5) 
MeV. At the peak, rjn, up and ao(980)er decay intensities are in the ratios 1 : 9.2 ± 0.8 ± 1.3 : 3.1 ± 0.2 ± 0.9. 
There is no evidence for a separate ao near 1300 MeV claimed by Obelix. Parameters of ao(980) are updated 
to M = 987.4 ± 1.0 ± 3.0 MeV, g 2 (yn) = 0.164 ± 0.007 ± 0.010 GeV 2 , g 2 {K K) / g 2 {t)-k) = 1.05 ± 0.07 ± 0.05. 
Its dominant second sheet pole in the KK channel is at (989 ± 1 ± 5) — i(40 ± 2 ± 4) MeV. Finally, the nature 
of the prominent J PC — h — > up signal in upn data is also clarified. 
PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 



1 Introduction 

The main objective of this work is to re-examine the parameters of ao(1450). It appears in the summary table of 
the Particle Data Group [1] although it has been observed decisively in only one set of data, from Crystal Barrel 
on pp — > i]ir°ir° at rest [2], [3], [4], [5]. There is further evidence from other data of the same experiment. It has 
been confirmed in the up channel in CL>7r + 7r~7r° data at rest [6]. Its KK decays are observed in (K1K ± )tt t data 
at rest [7]. There is also evidence for it in pp — > (ryV )^ [8] and pp — > (?77r + 7r _ 7r ± )7r T [9] at rest. 

Its branching ratio to up in Ref. [6] is a factor ~ 11 larger than to rjir (and will be revised here slightly). The 
fact that its branching ratio to t/tt is < 10% explains why ao(1450) has been elusive in data for n~p — > mrn. An 
important point is that the phase space for the up channel has a rapid s-dependence, so it is inappropriate to 
fit the ao (1450) with a Breit-Wigner amplitude of constant width, as was done in the early work. Its line-shape 
and the relation between magnitude and phase are affected strongly by dispersive effects, which are treated fully 
here. Attention to this detail improves considerably the stability of fitted parameters and makes the signal in 
both rjir and up much clearer. 

The ao (980) is examined along the same lines. The dispersive effect due to the opening of the KK channel 
plays a critical role and again improves the quality of the fit to r)ir°ir° data. Adler zeros are included into its tjtt 
and rf it decay channels. For these reasons, its parameters change significantly from earlier work. 

The dispersive effects may be unfamiliar to experimentalists, though well known to theorists since the 1950's. 
Experimental analyses have conventionally been done with a Breit-Wigner amplitude with denominator 

D{s) = M 2 -s-iY,9jPM- (1) 

3 

Here Pj(s) is the phase space for each decay channel j as a function of invariant mass squared s, possibly including 
a form factor. The gj are coupling constants to each decay amplitude. Let us write 

D( S ) = M 2 - S -^n i ( S ) (2) 

3 

with ImII(s) = g?pj(s). Because scattering amplitudes are analytic functions, any s-dependence of ImII(s) 
necessarily leads to a term in D(s) given by 

M*.^-/"^. ( 3) 

n Js(thr) S - S 

Here P denotes the Principal Value Integral and Sthr is the value of s at threshold. This is known as a dispersive 
contribution. It is equivalent to evaluating loop diagrams. If p(s) changes rapidly, as it does at the opening of 
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a sharp threshold, the dispersive term becomes dominant and affects the parameters of the resonance strongly. 
Fig. 1 below illustrates the result for ao(980). There is a prominent cusp in RcIIkk(s), centred at the threshold. 
It plays a major role in locking the resonance to this threshold [10]. One objective of the present work is to 
rchne the parameters of a (980) to include this effect. 

Consider next a (1450). The tup threshold is quite sharp and has a large effect on the line-shape near 1450 
MeV. There is a cusp at the cup threshold which also acts as an attractor. This may be the reason that ao(1450) 
is higher in mass than / (1370) and ifo(1430). 

Section 2 reviews dispersive effects. In principle they apply to all resonances. Fortunately, resonances with 
broad thresholds may be approximated by the pole term alone and this will be demonstrated here for a2(1320). 
There may be small residual effects far from resonance, but in practice these effects are tolerable at present. 

Section 3 discusses fits to ?77r 7r and W7r + 7r~7r° data, hence parameters of a (1450). In the present work, 
the widths of ao(1450) to KK and rjir are small, so there is no longer significant overlap between these two 
resonances and therefore little correlation between their parameters. An incidental feature of the re-analysis of 
cj7T + 7r~7r° data is an improved understanding of the large J PC = h ujp signal observed there. 

Section 4 gives results for ao(980) and Section 5 discusses KK coupling of ao(1450) and parameters of ao(980). 
The data on K ( [K ± 'k t do not give an accurate determination of their coupling to KK, but agree within sizable 
errors with the better determination from rjTr°TT° data. Section 6 summarises conclusions and makes some remarks 
on further desirable work. 

2 Technicalities of the dispersive terms 

As an introduction, let us consider ao(980) — ► KK. Mass differences between K + K~ , K° L K^ and K°K° will be 
ignored here because their separations are smaller than mass resolution in data to be fitted. There is a further 
reason. The VES group has very recently presented data showing that the /i(1285) decays to 3tt [11]. This 
violates isospin conservation and may well arise from mixing between do (980) and /o(980) due to mass differences 
in the KK thresholds. Consideration of this problem requires a combined analysis with data on / (980). It is 
necessary to take one step at a time and defer this for the present, though one should bear in mind there may 
be some small effect on parameters fitted to ao(980). 

Ignoring mass differences, pkk = yl — im^Js, where ttlk is the mean kaon mass, 495.663 MeV. As s — > oo, 
the phase space factor — > 1. Without any form factor, the dispersion integral of Eq. (3) diverges. Therefore a 
form factor 

F KK = exp(-afc 2 ) (4) 

is used to multiply Qkk ■ Here k is the momentum of each kaon in the KK rest frame. This well known form 
factor assumes a Gaussian source with RMS radius R given by a — R 2 /6. It turns out that the same value of 
a succeeds in fitting all resonances and avoids a multiplicity of parameters. It optimises at a = (2.0 ± 0.25) 
(GeV/c)" 2 , corresponding to R = 0.68 ± 0.04 fm. 




Figure 1: KK phase space p'(s), normalised to 1 at its peak and RcII(s) (full curves). Dashed curves show 
results for a = 2.5 (GcV/c)~ 2 . 

Library subroutines are available for evaluating the Principal Value Integral. Results are shown in Fig. 1 for 
the product p' — Pkk{s)F^ k (s) and for RcII(s) . For display purposes, the normalisation is chosen so that p' 
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including the form factor peaks at 1. There is a large cusp in the real part, somewhat larger than the peak of 
p' . The dashed curves show results with a = 2.5 (GeV/c)~ 2 . In the mass range where the ao(980) is strong, the 
sensitivity to a is quite small; it comes into play only in the wings of the resonance, where other resonances may 
mask its effects. 

There are two practical points concerning Rell(s). Although it is responsible for attracting do (980) to the 
KK threshold, it is convenient to make a subtraction in the Breit-Wigner denominator on resonance: 



Secondly, it is convenient to evaluate the dispersion integral as a 2-dimensional array against s and a. A simple 
sub-routine interpolates in this table. When fitting data, one can then optimise a, M 2 and g 2 with a standard 
optimisation program such as Minuit. The dispersive term is proportional to g 2 . Just below the KK threshold, 
it varies as g|fi<-(4m|- ~ s )/ s ( see algebra in Ref. [10], Eq. 15). This resembles the term (M 2 — s) in D(s). 
Consequently M and g 2 become strongly correlated unless there are data determining g 2 separately for every 
channel. 

Because of the correlations, the convergence of the fit it rather poorer than for a simple Breit-Wigner resonance 
of constant width but still adequate. It is in fact better to let the programme optimise the parameters. The 
alternative, a grid search over M and g 2 , is subject to the correlations between them. Standard optimisation 
programmes work with eigenvectors and circumvent the correlations. 

A general point is that all resonances are subject to opening thresholds, hence dispersive effects. However, it 
fortunately turns out that for broad thresholds the net effect of the dispersive terms becomes small within one 
full width of the pole. The data can then be parametrised directly in terms of the pole term oc l/(s — s po i e ). 
This conclusion emerged from work on /o(I370) concerning aa and pp thresholds [12]. For broad thresholds, 
the dispersive terms have significant effects only far from the pole, and may not be trustworthy there because 
of uncertainties in form factors. In the present work, thresholds for ai(1260) — > (pit)l=i, 02 (1320) — ► (pit)l=2, 
{t)tt)l=2 and (KK)l=2 have been treated fully using the dispersive term. These thresholds open fairly gently 
because of the centrifugal barriers for orbital angular momentum L in the decays. The conclusion is the same 
as in [12]: the line-shapes of these resonances are affected rather little except for their tails. 

However, in fitting upir data, the rather sharp top threshold does affect the fitted resonances quite strongly. 
The cusp at the up threshold is broadened by the line-shape of the p. This line-shape may be included in the 
evaluation of top phase space, then the dispersive effect can be evaluated from the phase space. Suppose as an 
example ao(1450) — > up, followed by p — > mr. The 3-body phase space for up is given by the integral 



where T is the Breit-Wigner amplitude for the p. Also s refers to the ao(1450) and k to the momentum of the u 
or p in the ao rest frame; s\ and k\ refer to the p and the momenta of the pions in its rest frame. When there 
is angular momentum in the decay to pu, a centrifugal barrier needs to be included. 

3 Fits to pp — > r]7T°7T° and utt + 7t^7t 

The ao(1450) was discovered in Crystal Barrel data for pp — > ?77r 7r at rest [2], [3], [4], [5]. It also appears in 
the up channel in upir data at rest [6], and in the KK channel in K^K^n^ data at rest [7]. The latter will be 
discussed in Section 5, but it turns out that the systematic error in its coupling to KK is rather large. Its KK 
coupling is consistent with the SU(3) prediction, and will be fixed to that value. The same applies to the rfir 
coupling. An analysis of data on pp — > ri'ir°ir° at rest gave results consistent with this prediction [8]. The effects 
of both KK and rj'n channels on the line-shape of a (1450) are similar to r/n and quite small. 

A preliminary comment is required on the fit to data for pp — > a>7r + 7r~7r , discussed in sub-section 3.9. The 
earlier publication did include dispersive effects. The fit to these data changes rather little in the combined fit 
with 777T 7r° data. The main improvement to parameters of a (1450) comes from the rjir a ir° data. 

The decays of consequence for the line-shape are rjn, up and ao(980)cr. If up is the only strong decay channel, 
it restricts the maximum possible full-width severely. This is because the rapid increase in p up with s inflates 
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the Breit-Wigner denominator and cuts off the upper side of the resonance. The consequence in Ref. [6] was 
that the optimum fitted full- width was only 110 MeV. 

This width is too small for good consistency with the r)ir°ir° data. Some other broad threshold is needed for 
an acceptable fit. This is provided by the ao(980)cx decay. Data for pp — > T/7r + 7r~ 7r~ tt~ at rest [9] were found 
to contain some ao(1450) signal in t]tt + t:^t: + , improving log likelihood by 32 for 2 fitted parameters; this is 
statistically > 7 standard deviations. However, there was no optimum when the mass and width of ao(1450) 
were scanned. The branching fraction for the rjAir final state is a factor 14 larger than for rjinr, with the result 
that the allowed branching fraction of ao(1450) — » ao(980)cr could be as much as 4.3 times that of t]tt. It now 
turns out that including the ao(980)<7 threshold supplies the required broad component in do (1450) decays and 
improves markedly the fit reported here. The required branching ratio to r)Air is only slightly smaller than that 
fitted in Ref. [9] , so it appears to be a genuine signal. 

Ideally the t]tt + tt~ tt + tt~ data should be included in the present fits. Unfortunately those data have been lost, 
so this is not possible without major work reprocessing them from raw data. This is not worth the effort, since 
they did not constrain the mass and width of ao(1450). All that matters is the magnitude of the fitted signal 
and the upper limit on the ao(980)cr branching fraction; these can be taken from the earlier publication. 




3.1 Features of the r]ir TT data 

The Dalitz plots for data in liquid and gaseous hydrogen are shown in binned form in Figs. 2 and 3. There 
are ~ 280,000 events in liquid hydrogen with experimental background < 1%. A minor detail is that any bins 
overlapping the edges of the Dalitz plots have been removed from Figs. 2 and 3 and the fits. There are also 
some further bins immediately adjoining edge bins and showing questionable behaviour. This can arise if an 
event lies outside the true Dalitz plot before the kinematic fit. That fit enforces the constraints of energy- 
momentum conservation and the masses of 7r°, r) and u>. It pulls events inside the Dalitz plot, but there is some 
tendency for them to congregate towards the edges. These bins are easily identified and removed because the fit 
is systematically lower than data. A total of 18 out of 3582 bins are removed for this reason, though effects on 
fitted parameters are tiny. 

Statistics for rjirir are so high (~ 280, 000 events) that it was not possible to equal those statistics in the Monte 
Carlo simulation. (Only a few per-cent of events survive the data selection). It is assumed that the acceptance 
is uniform, in accord with observations for present data, pp — > 777771"° and 37r°. The final fit has a x 2 of 2.9 per 
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bin. A similar value was obtained in fitting pp — > 3ir° [12]. Examination of the present fit reveals no systematic 
deviation across the Dalitz plot associated with fitted components. There are possible slow variations with \ 2 
up to 10 which could be associated with small systematic effects in the slowly varying 7T7t S-wave or alternatively 
could arise from small variations in experimental acceptance. Any departure from uniform acceptance over the 
width of ao(1450) has an effect much smaller than errors. However, it has been necessary to scale statistical 
errors to account for the mean \ 2 per bin. 
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Figure 3: The Dalitz plot for pp — > t/t: ^ at rest in gaseous hydrogen. 



Figs. 2 and 3 resemble one another closely, showing that the effects of P-state annihilation are small. Both 
sets of data are fitted fully, and the final fit gives 7.4% P-state contribution, in close agreement with an earlier 
determination [13]. This come mostly from 3 Pi and 3 P 2 -> a 2 (1320)7r and a 2 (1700)7r, and 3 Pi rja and 77/0 (980). 
The former plays an important role in fitting the angular dependence of the prominent a 2 (1320) bands. The 
effect of the latter two components is visible along the /o(980) band, where interferences between /o(980) and 
a affect the apparent width of the /o(980) in the data. Ultimately P-state annihilation has little effect on fitted 
parameters of either ao(980) or ao(1450). P-state production of ao(1450) is inhibited by a centrifugal barrier and 
makes only a very weak contribution (0.27%). 

The two a 2 (1320) bands interfere constructively at the upper right-hand edge of the Dalitz plots. Interfer- 
ence between the two a 2 's builds a bridge between them along this edge. The bands appear to be not quite 
vertical/horizontal. In the analyses of the 1990's, this deviation was fitted by a broad ?77r P-wave resonance with 
ill-defined mass and a large width of ~ 600 MeV. Those parameters are inconsistent with what is now known 
about the ?77r P-wave. The current fits are made with the tti(1400) parameters fitted to Crystal Barrel data on 
pn — > TT~TT°ri [14]. In those data, there is a significant P-state contribution because the process 1 P\ — > 7T7Ti(1400) 
goes via the S-wave. In present ?77r 7r data, there is now a small (0.9%) P-state contribution from 7r(1400) and 
only 0.6% in S-state annihilation. 

Another distinctive feature of the Dalitz plots is a sharp 'edge' in r/7r coinciding accurately with the KK 
threshold. This is due to the the opening of the KK threshold for do (980). At this threshold, pkk changes from 
real to imaginary as one crosses the threshold from above to below. The amplitude for ao(980) therefore turns 
in phase by 90°. Consequently, interference with the rnr S-wave changes dramatically. The precise form of the 
'edge' is therefore sensitive to the relative coupling of ao(980) between r]ir and KK. 

Fig. 4 shows mass projections for 7r°7r° and r]ir° in liquid hydrogen and the fit. Fig. 4(b) is the easier to 
understand. The first (left-hand) peak is a reflection of the a 2 (1320) at the left-hand side of the Dalitz plot. 
The sharp rise to the second peak is caused by the 'edge' due to ao(980) and its interferences with a 2 (1320). 
The third peak is directly due to a 2 (1320). The sudden drop at high mass is due to the a 2 (980) 'edge' crossing 
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Figure 4: Mass projections for (a) tttt and (b) rjn for pp — > t]tt tt at rest. Points with errors show the data; 
histograms show the fit. 



the right-hand side of the Dalitz plot. The quality of the data (and fit) illustrate the information available on 
ao(980) and its coupling to KK. Note that the ao(1450) is not directly visible in Fig. 4(b). 

In Fig. 4(a), there is one high point at = 0.44 GeV 2 . It does not correlate with anything and appears to 
be a statistical storm. The first peak to its right is again due to ao(980) and its interferences with ct2(1320) and 
the tttt S-wave. The second peak is a reflection of 02 (1320) on the lower side of the Dalitz plot. 

The cusp in ao(980) at the KK threshold is sufficiently narrow that it is necessary to fold in the mass 
resolution for bins adjoining the KK threshold. The mass resolution is a Gaussian with a a of 9.5 MeV. This 
number is derived from data on pp — > tt 77' (958) [15], where the fitted width of the r\' is readily measured. The 
folding is done using Gaussian 12 point integration over the bins concerned. 

3.2 The ttty S-wave amplitude 

This is the third major component in the rjir°ir° data. Since the earliest publications in the 1990's, our knowledge 
of the tttt S-wave amplitude has improved greatly. Today, the a pole is well known from (a) the BES2 data on 
J/'fy — > ujir + Tr~ , where it produces a strong peak at ~ 500 MeV [16], (b) the calculations of Caprini, Colangelo 
and Leutwyler using the Roy equations to constrain the s-dependence of the elastic amplitude [17]. 
The tttt elastic scattering amplitude may be written in the form 

f el (s)=N(s)/D(s), (7) 

where N(s) is real and must be equal to — ImD(s) below the KK threshold. In a production reaction, D(s) 
must be the same for the a pole as in elastic scattering (Watson's theorem [18]). However, N(s) is allowed to be 
quite different between production and elastic scattering [19]. The strong peak close to 500 MeV in BES data for 
J /IS? — > ujtt + tt~ is fitted accurately taking N(s) to be constant. There is then accurate agreement [20] between 
the pole observed in these data and the elastic phase shifts predicted by Caprini et al. 

The recent fits to data on pp — > 3tt° [12] require a 2-component form for the S-wave production amplitude: 

fprodn = Ai/ ei (s) + A 2 /D(s), (8) 

where Aj.,2 are complex coupling constants: i.e. a coherent sum of the elastic amplitude and the pole term. This 
2-component prescription also fits the r)ir°Tr° data, with different A to those for pp — ► 3tt°. This prescription will 
play an essential role throughout the present work, including the fit to data on pp —> K\K^tt^ where both the 
k pole and the clastic Kit amplitude contribute. Similar variations of the tttt S-wave amplitude are well known 
in decays of T' to Ttttt and are discussed in detail by Simonov and Veselov [21]. 

Two alternative prescriptions are available for the tttt S-wave, from Refs. [20] and [12]. The latter is fitted to 
data for pp — ► 3tt°, where the tttt mass range extends to 1.74 GeV. For pp — > 7771"° 7r°, the tttt mass range stops at 
1.329 GeV. The Att inelasticity is quite small up to this mass. The two alternatives lead to only minor differences 
in the quality of fit to pp — > r)TT°TT°. The first prescription is simpler and faster and is used for final fits. 
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A question arises whether to assume the a (980)er channel is produced via the a pole or the elastic tttt 
amplitude. If the latter is used, a (980)<7 phase space rises too slowly to have much effect over the mass range 
of ao(1450). Some production of the tttt S-wave via its pole term is needed and is what is used here. It is also 
what was fitted to r]4:TT data. The mean mass of the a is then ~ 470 MeV and the full width is ~ 500 MeV. The 
a <7 phase space is then substantial at 1450 MeV. 




m (GeV) m (GeV) 

Figure 5: (a) Phase space for top (full curve) and ao(980)a (dashed); (b) contributions to ReII(s) for a (1450). 

Fig. 5(a) shows cop and ao(980)a phase space. They both peak in the mass range 1.65-2.0 GeV because 
of the form factor. Fig. 5(b) shows the subtracted form m(s) = ReII(s) — ReII(M 2 ) for ao(1450) with the 
normalisation of final fits. To a first approximation, they are proportional to the gradient of phase space. Their 
sum is roughly half the magnitude of M 2 — s, so their effects on the line-shape of the resonance are quite large. 

3.3 Treatment of / (980) 

A further element required to fit the r]Tr n data is from / (980). A full re-analysis of /o(980) parameters 
including the dispersive effect is a major undertaking requiring fits to the many sets of data in which it is 
prominent. For present data, the line-shape of /o(980) is not critical. Fig. 6(b) below will illustrate the blurring 
of the threshold cusp in a (980) due to mass resolution. The blurring is even more severe for /o(980), which has 
a full- width at half-maximum of only 34 ± 8 MeV [22], compared with the mass resolution of ±9.5 MeV. The 
mass resolution is folded with the line-shape of /o(980) in the fit to data. The effect of the threshold cusp is a 
marginal decrease in the width compared with the BES parametrisation. This small perturbation has negligible 
effect on the parameters fitted to ao(980) and a (1450), because information on these two resonances comes from 
regions of the Dalitz plot having only modest overlap with / (980). 

3.4 Fits to a (1450) 

Fits have been made simultaneously to the four set of data listed in the Abstract. The precise formula fitted 
to a (1450) needs discussion. A form factor is needed in calculating dispersive terms, in order to make the 
dispersion integrals converge. However, as Fig. 5(a) shows, the form factor plays a strong role only above 1650 
MeV, well above the ao(1450). The form factor is therefore an unnecessary elaboration over the mass range 
covered by ao(1450). For simplicity, it is therefore dropped in the amplitude fitted to data. 
The ao(1450) amplitude may be written 

/(1450) = 1/[M 2 - 8 (U(s) U(M 2 )) - i£>fo(*)], (9) 

3 

where the sum runs over 7777, KK, rj'ir, top and a (980)cr channels. A slightly rearranged formula will be given 
later in the light of observed results. Values of g 2 for KK and r]'ir at the peak (i.e. near the pole) will be 
fixed to SU(3) predictions, which depend on the angle <j> — 54.7° — 6ps, where Ops is the pseudoscalar mixing 
angle. Values of <j> may be obtained from analysis of radiative decays of vector (V) and pseudoscalar mesons 



7 



(P) mesons. Escribano and Nadal analyse all existing data and conclude there is no significant evidence for a 
gluonic component in 77 or rj [23]. It seems prudent to use results without that component. They then find 
4> = (41.5 ± 1.2)°. Thomas does a similar analysis with an identical conclusion [24]. Data on J/^> — ^ VP decays 
also give a less precise result: <p = (40.5 ± 2.4)° [25]. The weighted mean </> = (41.3 ± 1.2)° will be used here. 



Then 

Sfyjg^ = tan 2 ^ = 0.772 ± 0.068, (10) 
9 2 kk/9% = 1/(2 cos 2 </>) = 0.886 ±0.034. (11) 

3.5 Treatment of branching ratios 

The relative value g^ p /g^ is obtained from relative branching fractions of ao(1450) in riirir and ujpir data. 
However there are two points which need to be taken into account. 

Firstly the observed branching fractions for each resonance in pp data must be obtained by folding the phase 
space factors for each channel with the line-shape of the resonance, using integrals of the form 



/ 



f dsg](s) Pj (s)k' 
3 ~ 1 \D{s)f ' 



The factor k' is the momentum with which the resonance is produced in pp — > -tt + ao(1450); it allows for the 
phase space corresponding to the length of the ao(1450) band as a function of s on the Dalitz plot. 

Secondly, there is an important point of principle concerning how to account for interferences. Data are 
fitted including all the interferences, not only between different resonances but also including, for example, two 
a2(1320) appearing in r)ir°ir° data. The coupling constants A are determined by the fit; but then, for use in Eq. 
(12), intensities of individual components must be evaluated from these A without the interferences. The two 
a 2 's contribute 30.3% of r]TT Tr data, but 5.8% of this arises from interference between the two bands. 

There are even larger effects for ao(1450). In rjTnr, there are constructive interferences between the two 
ao (1450). Including interferences, they contribute 5.44% of the cross section, but without interferences, this drops 
to A = 3.48%. In cj7r + 7r~7r° data, there are three charge states for ao in the amplitude (a^ir^ — aQ7r + Oq 7r + ), 
where signs arise from isospin Clebsch-Gordon coefficients. There are therefore some destructive interferences. 
With this interference included, the a 's contribute 3.49% of cj/ot data, but without them B = 4.86%. The ratio 
B/A determines Sg^pI^p/g^I^. For the ao(980), interference effects are quite small, because the peak of the 
resonance is narrow. 

It is necessary to arrange, iteratively, that the fitted branching ratio between rjir and wp signals is consistent 
with the fitted value of a 2 /a 2 . Table 1 lists the percentages of the signals fitted to i]tt°tt data including 
interferences. These do not add up to 100% because of interferences. 



signal 


Percentage 


7T7T S-wave 


9.4 ±0.3 


/o(980) 


11.7±0.2 


a (980) 


12.8 ±0.2 


a (1450) 


5.44 ±0.25 


a 2 (1320) 


30.3 ±0.2 


a 2 (1700) 


6.2 ±0.2 


7^ (1405) 


0.6 ±0.2 


P-states 


7.4 ±0.5 



Table 1: Intensities of signals fitted to pp — ► ?77r u 7r" data. 



3.6 Comments on the fit to a (1450) 

The fit to a (1450) improves significantly compared with work in the 1990's where a Breit-Wigner resonance of 
constant width was assumed. Parameters of a (980) and a (1450) were correlated significantly in that early work, 
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because the large and constant width of a (1450) made it overlap a (980) significantly. The mass of a (1450) 
could move between 1450 and 1510 MeV as T was varied. Now the width of ao(1450) near 1 GeV is restricted to 
rjir and KK and the width to r/ir is only 19 MeV at a mass of 1 GeV. The result is that ao(980) and ao(1450) are 
now almost uncorrelated. With the s-dependent forms used here, the peak position is very stable in the range 
1440-1460 MeV, with an optimum at 1448 MeV. 

Furthermore, the fitted ao(1450) — > rjir signal increases from 3.0% to 5.4%. The data clearly prefer the 
s-dependent form. The full width of ao(1450) at half maximum decreases substantially from the 265 ± 13 MeV 
quoted by the Particle Data Group [1] to 192 ± 9(stat) ± 9(syst) MeV. This is inevitable in view of the rapidly 
increasing oo(980)(j and Lup signals, which make the Breit-Wigner denominator cut off the line-shape at high 
mass. The PDG value is subject to serious systematic error from the assumption of constant width. 




mass (GeV) mass (GeV) 

Figure 6: (a) Natural line-shape of a (1450) (full curve) and with the effect of pp phase space included (dashed 
curve); the dotted curve shows a Breit-Wigner resonance of constant width agreeing at half- height with the full 
curve, (b) Natural line-shape of do (980) (full curve) and including the mass resolution of the Crystal Barrel 
detector (dotted). 

The ao(1450) line-shape for an isolated resonance, i.e. without the factor k! of Eq. (12), is shown by the 
full curve of Fig. 6(a). What is plotted is \1/D(s)\ 2 , i.e. ignoring any phase space effects in the numerator of 
the amplitude. The line-shape observed in pp data including the factor k! is shown by the dashed curve. A 
Breit-Wigner line-shape with constant width is shown by the dotted curve, agreeing at half-height with the full 
curve. The true line-shape is asymmetric because of the rising phase space for uip and ao(980)cr and also because 
of the dispersive term IT(s) in the Breit-Wigner denominator. 

Fig. 6(b) shows the line-shape of ao(980) without and with the effect of mass resolution of the Crystal Barrel 
detector. The ao(980) is cut almost in half by the opening of the KK threshold. At this threshold, the line-shape 
drops rapidly because of the KK width in the Breit-Wigner denominator. Many theorists base calculations on 
the 50-100 MeV width of the ao(980) quoted by the PDG. This is the full width at half-maximum. The values 
of g 2 and g 2 KK are both ~ 160 MeV, comparable with other resonances. 

The Argand diagram for the coherent sum of ao(1450) and ao(980) is shown in Fig. 7, excluding the effect of 
mass resolution. The maximum amplitude for ao(1450) is at 1448 MeV, where the phase of the ry7r amplitude is 
only ~ 50°. The phase goes through 90° only at 1536 MeV. This is the mass M in the Breit-Wigner denominator. 
The full curve of Fig. 8 shows the Argand diagram of ao(1450) drawn from 1/ D{s) alone. It appears to lie on 
its side because it is cut off at high mass by the rapid increase of agcr and up phase space. The value of M in 
the Breit-Wigner denominator is a derived quantity, rather strongly dependent on g 2 (u>p) and <7 2 (ao(980)er) and 
their form factors. Accordingly, M has a factor 2 larger error than the peak mass, which responds directly to 
the magnitude of the rjir amplitude. 

An important check is whether ao(1450) really requires a resonant loop like that shown on Fig. 7. The first 
check is to replace the resonant form with its absolute magnitude, deleting its phase variation. As expected, x 2 
increases by 297.3 (after renormalising to allow for the fact that x 2 is 2.9 per data point); this is a 17 standard 
deviation effect. 
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Figure 7: Argand diagram for the coherent sum of ao(980) and do (1450), excluding the effect of experimental 
mass resolution. Masses are shown in GeV. 

A more delicate check is to break the mass range from 1315 to 1675 MeV into 30 MeV bins and optimise the 
ao(1450) signal in each bin. The result is compared with the Argand diagram of ao(1450) alone on Fig. 8. The 
individual bins follow the expected loop closely up to 1540 MeV; dotted lines show the movement of each bin 
from the analytic formula. Up to this mass, there is large interference between 02(1320) and ao(1450), providing 
strong constraints on its phase variation with mass. Above this, three of the four remaining points show a rather 
large scatter. Above 1560, the effect of 02(1700) becomes more important than that of d2(1320). Final fits use 
Crystal Barrel parameters for a 2 (1700): M = 1660 MeV, T = 280 MeV [26]. The main problem is that the 
corner of the Dalitz plot above m,„ = 1560 MeV is a cramped area in which to separate spin zero components 
from spin 2. The a2(1700) has significant contributions from all three initial states So, 3 Pi and 3 P2- These 
allow it to simulate a spin contribution to some extent, despite the existence of the data in gas which help 
determine P-state contributions. The most likely explanation of the discrepancies above 1560 MeV is a poor 
separation between a 2 (1700) and do(1450). 
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Figure 8: The Argand loop of ao(1450) (full curve) compared with magnitude and phase for individual bins 30 
MeV wide, centred at 1.33 to 1.66 GeV. Dotted lines show how individual bins move from the analytic formula. 

Finding the pole position of ao(1450) requires parametrisations of a^a and pu phase space and the dispersive 
term in the Breit-Wigner denominator. This has been done with three alternative parametrisations for each of 
the three terms. Formulae are chosen with good convergence properties for complex s, i.e. powers of s confined 
to the denominators of the formulae. All combinations of the formulae agree within ~ 3 MeV for both real and 
imaginary parts of the pole, showing that systematic errors for the extrapolation are well under control. The 
pole position is 1432 — i98 MeV; the main systematic errors arise from the mean mass and width of the peak. 

Fig. 9 shows as dashed curves the intensities of ao(1450) — > ao(980)er and uop as they appear in production 
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from pp. The a (980)cr decay peaks at 1458 McV and the ujp decay peaks at 1476 MeV. Curves are normalised 
to 1 at their peaks. Full curves show the peaks for an isolated resonance without the limitation of pp phase space 
for production. The ao(980)a peak is then at 1467 MeV and the top peak at 1485 MeV. 



a (980)a 




b) up 




mass (GeV) mass (GeV) 

Figure 9: Line-shapes of ao(1450) as it appears in decays to a (980)cr and cop, all normalised to 1 at the peaks. 
Dashed curves include the phase space for production in pp annihilation; full curves are for an isolated resonance 
without the limitation of the production process. 



3.7 Fitted parameters 

Table 2 collects results for ao(1450) from the final fit. At this point, it is necessary to present a more convenient 
formula for ao (1450) than Eq. (9), and the rationale behind it. 

The basic points spring from the fact that there is a pole at 1432 — i98 MeV. If one knew in advance how 
parameters vary between the pole and the physical region, it would be best to write the formula directly in terms 
of the pole and its residues, which express its coupling to every channel. That is not the case, so the closest 
approach is to write the formula in terms of the nearby peak mass, m p = 1448 MeV and widths to each channel 
at this mass, together with their s-dependence: 

/(1450) = 1/[M 2 - s-i[U(s)-n(M 2 )]-im p J2 T j( s )h (13) 

3 

mpTjis) = g^Pj(s). (14) 

This form is close to that for a Breit-Wigner resonance of constant width and is closely related to observed 
branching ratios between channels. 

The branching fractions for an isolated resonance are given by integrals of the form 

f 9jPj(s) ds 
J \D( S )\ 2 ■ 

If pj, hence Fj(s), were to vary linearly with s, the variation of branching fractions would cancel between upper 
and lower halves of the peak. It turns out that this cancellation works fairly well. This form of parametrisation 
gives a clear insight into the way the fit responds to each parameter. 

However, one important point emerges. The sum of the widths at the peak comes to 345 MeV, considerably 
larger than the observed full width of the peak, 192 MeV. The reasons for this are straightforward. On the lower 
side of the peak, r(ao(980)(r) and T(ujp) are small, and the amplitude falls rapidly because the remaining width 
to ?77r, KK and r]'ir is small. On the upper side of the peak, r(a (980)cr) and T(cop) rapidly become large and 
dominate the denominator, cutting off the line-shape l/\D(s)\ 2 quickly. The line-shape of ao(980) in Fig. 6(b) 
serves as a second example. The upper part of the peak is attenuated rapidly by Tkk(s)- The lower part is not 
far from a Breit-Wigner resonance of constant width. 

The branching ratio r(ao(980)cx)/r(777r) is obtained from the combined fit. It is consistent with the magnitude 
of the a (1450) signal fitted to rj-K + n~ 7r+ tt~ data in Ref. [9] and with the upper limit of 4.3 established there for 
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1432 ± 13 ± 25 - i(98 ± 5 ± 5) 



(3.49 ±0.14 ±0.30)% 
(4.86 ±0.19 ±0.42)% 



(5.44 ±0.15 ±0.88)% 
(3.48 ±0.13 ±0.58)% 



23.7 ±0.5 ±2.0 
17.7 ±0.3 ±2.0 
11.4 ±0.2 ±1.5 
219 ± 18 ±24 



1448 ± 13 ± 25 
1536 ±20 ±30 
1424 ± 13 ±25 



2.3 ±0.2 ±0.6 
7.6 ± 0.6 ± 1.2 



73 ± 5 ± 20 



192±9±9 



Table 2: Results for a (1450) in units of MeV. The first errors are statistical and the second systematic. 



this ratio. The Table uses branching ratios to cop from the combined fit to cj7r + 7r _ 7r° data, discussed in detail 
in subsection 3.9. The Table quotes in lines 12 and 13 the ratio of branching fractions B R(u) p) / B R(r]ir) and 
B i?(ao(980)er) / B R(t]tt) as they appear in pp annihilation. These values are better defined than those for an 
isolated resonance because of uncertainty about its high mass tail, see Fig. 9. Note that B R(cu p) / B R(i]ir) is 
7.6 in Table 2, rather smaller than the ratio T(ujp)/T(r]Tr) = 9.2 at the peak of ao(1450). This is because the 
up signal in pp — > ao(1450)7r is inhibited at high mass by the available phase space, as illustrated in Fig. 9(b) 
by the dashed curve. For an isolated resonance, the top branching fraction is close to 9.2, but with an unknown 
error depending on form factors. 

In Table 2, the first errors are statistical and the second systematic. Strong contributions to systematic 
errors arise from uncertainties in branching fractions for pp — > ?77r 7r (±5.0%) and u>7r + 7r~7r° (±7.9%). However, 
the largest error arises from the fact that interferences within one set of data lead to a branching fraction 
a A*/i| 2 ; rather than J2 i \kifi\ 2 . Here, the sum is over resonances, A, are coupling constants and /j are 
amplitudes for each resonance. The second of these quantities is derived from A parameters fitted to the first, 
as explained above in subsection 3.5. 

There is a potentially large error from the interference between the two components making up the tttt S-wave: 
the a pole term and the elastic component. Fortunately, the r]TT ir data determine both relative magnitudes 
and phases of these two contributions quite well. However, it is necessary to add a systematic error to cover the 
change to the fit if a further term is added to the parametrisation of the tttt S-wave. Here, it is chosen to be the 
clastic amplitude multiplied by s. There is a further small contribution to systematic errors from perturbations 
when small components are dropped from the fits, e.g. the weak 7Ti(1400) contributions in both 1 So and 3 Pi 
annihilation. Finally, in view of the scatter of the last 4 points of Fig. 8 above 1540 MeV, a systematic error is 
included from changes in the fit if do (1450) is fitted to i]tt°it data only up to 1540 MeV. 

Table 2 includes systematic errors in fitting lutt + tt^tt° data. The evaluation of systematic errors for these 
data follows the same procedure as for 7/7r 7r°. The final systematic errors are added in quadrature. It is not 
correct to add them linearly, as is sometimes done. The derivation of the Gaussian error distribution depends 
on the convolution of many box-shaped distributions. 
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3.8 A disagreement with Obelix 

The Obelix group has published two claims to observe an ao decaying to f]n in the mass range 1290-1313 McV 
[27] [28]. Such a resonance should be very conspicuous in Crystal Barrel data through distinctive interference 
with a2(1320). The fits reported here have been repeated (i) using an ao in this mass range without ao(1450) 
and (ii) together with ao(1450). When ao(1450) is removed from the fit, x 2 (scaled to allow for the mean x 2 of 

2.9 per bin) is worse by 528 for a reduction of six fitting parameters. This is an 18 standard deviation signal. 
If its mass and width are moved down to the mass range 1200-1340 MeV with a width < 120 MeV, there is no 
optimum. Instead the fit moves in a few iterations towards parameters of ao(1450), whatever line-shape is used 
for a (1450). If an extra a is added in the mass range 1280-1340 McV, there is only a small improvement in 
X 2 and again no optimum for parameters in the range claimed by Obelix. The narrow width they claim <~ 80 
MeV is similar to that of a 2 (1320). It appears likely that their signal was confused with P-state annihilation to 
a2(1320). The P-state annihilation is precisely identified in present work from data in hydrogen gas. 

3.9 The fit to ujtt+tv "tt data 

There are 35,280 reconstructed events for these data with 8.4% experimental background, arising in the selection 
of the narrow to. The earlier analysis of these data is reported in detail in Ref. [6]. Dispersive effects were 
included fully and the new fit changes rather little. Table 3 lists the components in the fit and their significance 
levels, measured by changes in log likelihood when each component is removed from the fit and all others are 
re-optimiscd. Values of x 2 are twice those for log likelihood for the large statistics available here. 



Initial states 


Channel 


Intensity (%) 


ALn L 


l S 


a (1450)7r, 


3.5 


90 




6i(1235)p, 61 -» ton 


13.2 


361 




7Tl(1600)7T, 7T1 — > [oi7r] L=0 


6.6 


71 




a 2 (1320)7r 


2.0 


18 




a 2 (1660)7r 


2.4 


36 




All 7T 7r 


16.9 


505 




ai(1260)7r 


0.8 


48 




oi(1260)w, ai — > pn 


23.9 


377 




ai(1640)7r 


5.1 


271 




ai(1640)7r, ai — > biir 


1.7 


45 




7ri(1600)7r, TTi -> [o l7 r] L=0 


2.5 


50 




w(1420)tt, w(1420) -» iva 


1.2 


8 




w(1420)tt, w(1420) -> bm 


3.3 


5 




61 (1235)o-, 61 -» L07T 


1.6 


57 




p(1450)a, p —> 6itt 


0.6 


94 




7T 7r 


9.5 


83 


6 Pi, 6 P2 


a 2 (1320) 


2.5 


38 




a 2 (1660) 


5.7 


45 



Table 3: Percentage contributions of each channel after the background subtraction. Decays are to uip unless 
stated otherwise. The final column shows changes in log likelihood when each channel is removed from the fit 
and remaining contributions are re-optimised. 

In the earlier work, there was a very marginal signal due to ai(1260) — ► uip, which improved log likelihood 
by 24. In the latest work, it improves log likelihood by only 6 and is omitted from the fit. Likewise the earlier 
work included a rather marginal signal for 7r(1600) — > (6i(1235)7t)l =2 . This contribution is now small and is set 
to zero. 

A further detail is that there are data for pp — > cj7r°7r 7r° [29]. The branching ratio for this channel is very 
small. These data constrain the magnitudes of the last four entries to Table 3 for 3 Si. Their phases are fitted 
freely. Two of them have only very small effects in the present fit. 
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Let us review the essential points of the analysis. There are three charge combinations of irir. As a result, 
individual resonances do not appear clearly in mass projections. It is necessary to rely on the amplitude analysis 
to locate magnitudes and phases from what it finds in 4-body phase space. That may appear questionable, but 
in practice works well. The fits to mass projections were shown in Fig. 2 of the earlier paper. The tiny changes 
in the new fit are hardly visible by eye and therefore the figures will not be repeated here. 

Secondly, angular distributions depend distinctively on spins. Consider ao(1450) as an example. The spin of 
the u) lies along the normal to its decay plane. The spin of the p is given by the vector (k\ — fc 2 ), where k\ t 2 are 
momenta of the pions from its decay. After Lorentz transformations to the top rest frame, the matrix element is 
given by the scalar product of these two vectors. This is highly distinctive. An elementary check on formulae is 
that all amplitudes are orthogonal. One can test how well quantum numbers are recognised by putting deliberate 
errors into formulae. Generally the result is that the amplitudes drop to small values. 

A third point is that the programme prints a matrix giving intensities of all components together with real 
and imaginary parts of all interferences. This identifies the important interferences. It is then easy to test the 
reliability of these interferences by plotting log likelihood against relative phases. 

The magnitude of the ao(1450) signal has decreased slightly from the earlier publication, but its significance 
level has improved. In the earlier work, log likelihood changed by 56 when a (1450) was omitted from the fit. 
Now it changes by 90. It is produced from the 1 S'o initial state where it interferes with a large and well identified 
6i(1235)p signal. 

There is a large signal in 3 5i annihilation from J PC = h . It peaks at 1480 MeV, quite close to ao(1450). 
One might worry that there will be cross-talk with ao(1450), despite the fact that the a (1450) is produced from 
the 1 Sq initial state while h is produced from 3 Si. There is no such problem. The h component may be 
removed completely from the fit without any significant effect on the fitted ao(1450) signal. In fact, the ao(1450) 
is insensitive to changes in all 3 Si, and P-state amplitudes. 
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Figure 10: (a) The Argand loop of ao(1450) in upir (full curve) compared with magnitude and phase for individual 
bins 60 MeV wide, centred at 1.345 to 1.645 GeV. Dotted lines show how individual bins move from the analytic 
formula, (b) The ratio of the ao(1450) amplitude in bins 30 MeV wide to that of the overall fit. 

A similar comparison has been made in Fig. 10(a) to that shown in Fig. 8. Real and imaginary parts of the 
a (1450) — > Lop signal are fitted freely in bins 60 MeV wide from 1.345 to 1.645 GeV. In Fig. 10(a), the amplitude 
is k' \jpup/ D(s); i.e. it allows for the phase space of the up final state and the phase space in the production 
reaction pp — > ao(1450)7r, proportional to the momentum k! of the ao(1450) in the pp centre of mass. There 
is no doubt that the data conform with a resonant circle, though errors are sizable. There appears to be some 
tendency for the data to require a larger amplitude than the overall fit. An alternative test is made by fixing the 
phase of the amplitude to that of the overall fit, but allowing the magnitude of the fitted signal to fit freely in 30 
MeV wide bins. Fixing the phase stabilises the fitted amplitude considerably. Results are shown in Fig. 10(b). 
There is now a reduced tendency for the fitted amplitude to be above the overall fit. The mean discrepancy is 
(13.7 ± 8.5)%. The systematic difference arises because the final fit is constrained to fit the line-width of the 
ao(1450). If the ao(1450) — > top signal is increased, it makes the line-width smaller; this was the problem with 
the first publication, Ref. [6]. The final fit is a compromise between fitting the line-shape and the magnitude of 
the ao(1450) — ► top signal. 
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In assessing the errors for T(up)/T(r]ir) in Table 2, the statistical error is taken from the 8.5% statistical error 
in the discrepancy of Fig. 10(b). This is quite close to the error derived from log likelihood in Table 3. The 
systematic error is derived from changes in r(wp)/r(777r) as the mass and width of ao(1450) are varied over the 
range of systematic errors in Table 2. 

One new point does emerge from a better understanding of dispersive effects. This concerns the large 
J p = h component. In the earlier work, attempts were made to fit it with 7r(1300) and a radial recurrence 
in the mass range 1600-1700 McV. However, the required signal for 7r(1300) was unreasonably large and would 
have required it to decay dominantly to up. Furthermore, the data still required a definite peak in the vicinity 
of 1500 MeV. A radial recurrenece so close to 7r(1300) would be surprising. 

The present work reveals a more sensible way of fitting the 0~ signal. The optimum fit is obtained with a 
broad resonance at 1540 MeV with Y = 590 MeV, plus a radial excitation of tt(1300) at 1732 ± 32 McV with 
T = 252 ± 30 MeV. The broad resonance is close to being a simple cusp at the up threshold. This solution is 
shown in Fig. 11. 
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Figure 11: (a) The line-shape of the J PC = 0~ + signal as it would appear in up without the limitation of 
production in pp — > upir but (i) including up phase space (full curve) and (ii) with this phase space factored 
out (chain curve) . The dotted curve shows the phase space for up including the limitation due to production in 
pp — > upir. (b) The Argand diagram; masses are shown in GeV. Units are arbitrary for display purposes. 



Fig. 11(a) shows as the full curve the intensity of the fitted up signal including up phase space, for a resonance 
'in free space', i.e. without the limitation imposed by production in pp — > upir. The phase space for pp — > upir 
including this limitation is shown by the dotted curve. What is actually observed in uir + ir~ir° data is then given 
by the product of the dotted curve and the full curve. The chain curve shows the result of dividing the full curve 
by up phase space. This result shows the line-shape arising from \Aa/ 'Da(s) + As / 'Db{s)\ 2 , where A and B 
refer to the two components and A are coupling constants. This is what one would see for an isolated resonance 
if there were no phase space factor in the numerator. The double peak near 1.45 GeV comes from the cusp + 
interference. The peak at high mass comes from the radial recurrence at 1732 MeV. 

Minor variants on this solution are possible because the up amplitude below 1450 MeV is small and poorly 
determined, particularly its phase. One should therefore not place any reliance on the threshold behaviour of 
the up amplitude below 1450 MeV. It is quite possible that the double peak at 1430 and 1480 is an artefact. 
However, there are two essential features which are unavoidable. The first is a peak at <~ 1550 MeV in up. This 
arises from the cusp at the up threshold. The second well determined feature is the phase advance > 180° from 
1510 to 1740 MeV. This arises largely from the radial recurrence at 1732 MeV. 

The present data are limited by the fact that production of 0~ from 3 Si is suppressed at the highest masses 
by the L — 1 centrifugal barrier for production. A quite significant signal is however observed also in P-state 
production with L = 0. This signal gives a reasonable determination of the mass and width of the upper 
resonance at 1732 MeV. 

There is some chance that this radial excitation corresponds to 7r(1800). However, using PDG parameters 
for 7r(1800), the fit is 2.8 standard deviations worse than with a free fit. The mass and width observed in the 
present fit correspond closely to those observed by Amelin ct al. in ir~ A — > uir~ir°A* [31]. The ideogram shown 
by the PDG for the mass of 7r(1800) has a double-humped structure. There is the possibility of a 0~ hybrid in 



15 



this mass range to accompany the 7r 1 (1600). So there is room for a conventional radial excitation around 1730 
MeV and a hybrid at higher mass. The 7r(1800) has decay modes suggestive of a hybrid. A fit using a mass and 
width from the higher lobe of the PDG's ideogram is worse than the free fit by 4.0 standard deviations. Further 
exploration of 0~ signals is needed in this mass range to resolve the current uncertainties. 




mass (GeV) Re f 

Figure 12: (a) The line-shape of a2(1320): full curve including s-dependence of the width and dispersive effects 
in full, the dotted curve showing a Breit-Wigner amplitude of constant width, (b) The Argand diagram. 

A final detail concerns 02(1320). Fig. 12 illustrates the small effect on a2(1320) of the s-dependence of the 
width and dispersive effects. The full s-dependence of decays to p7r, tjtt, KK and up is included. The full curve 
on Fig. 12(a) shows the fitted line- width and the dotted curve the line-shape of a Breit-Wigner amplitude of 
constant width agreeing at half-height with the full curve. There is little difference between them, showing that 
the slow s-dependence of the dominant pir channel has little effect. Fig. 12(b) shows that the Argand diagram 
follows a circle closely. 

4 Results for a (980) 

In the 1994 work [3], the a (980) was fitted with a Flatte formula with M = 999 ± 5 MeV, g 2 (r]Tr) = 221 ± 20 
MeV and g 2 (K K) / g 2 (i]tt) = 1.16 ± 0.18. The r/'ir channel was not included. These parameters now change 
beyond their errors because (a) the large dispersive cusp at the KK threshold is included, (b) Adler zeros are 
included in both rjir and rj'ir channels. These are at s = sa = m 2 — 0.5m 2 = 0.2905 GeV 2 for the rjir channel 
and at s = sa> = m 2 , — 0.5m 2 = 0.9078 GeV 2 for the r/'ir channel. Note that the latter is quite close to the 
resonance mass. 

The formulae used for coupling to rjir and rj'ir are 

9lM = g 2 m {Am 2 K ) S ; SA cxp(-2 a fc 2 ), (15) 

gl'As) = ^("V + m 7r ) 2 - V2 exp(-2a/c 2 ,), (16) 

ym,^' -\- m,,- ) sa' 

where k v ' is the momentum in the rj'ir rest frame. The value of g 2 , n is normalised at the rj'ir threshold. 
Below the thresholds for these processes, the Flatte formula has the sub-threshold analytic continuation 

9$>*Pv'n i9v>ir2\tf>\/s. (17) 

However, g 2 , v becomes real again for a mass < 0.83 GeV, due to the opening of the w-channel. It makes little 
sense to allow for this without including the dynamics of the w-channel process. Therefore g 2 ,^ is set to zero 
below s — sa>- The effect of the rj' channel on x 2 of the fit is quite small, except above the rj'ir threshold. Its 
coupling constant is fixed to that of the rjn channel by Eq. (10). Note that the assumption is made in Eq. (16) 
that the form factor docs not affect the ratio g 2 , v /g 2 ^ between these two thresholds. Because of the small effect 
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of the r/'ir channel on present fits, this assumption has little effect on \ 2 - However, it could matter if and when 
data become available directly on the rj'n channel. A detail is that the opening of the tj'tt channel is visible on 
the Argand diagram of Fig. 7, just below 1.1 GeV. 

For the KK channel, the Adler zero is far away at s = Q.hm 2 K and experience with the a amplitude is that 
the factor (s — SA)/{^'m 2 K — s) needs to be multiplied by an exponential form factor exp(— as), which prevents 
the amplitude rising indefinitely with s. Below the KK threshold, the KK channel has only an indirect effect 
on data in the r/n channel. Tests have been made with a variety of form factors. Within errors, the best fit is 
obtained with g 2 KK — constant below the KK threshold and this simple prescription has been adopted. Above 
the KK threshold, the factor due to the Adler zero is dropped and the form factor exp(— 2k 2 KK ) is used, with 
kpcK the momentum in the KK channel in GeV/c. 

4.1 Fits to 7]n°n° data 

The dispersive cusp locks the mass M of the amplitude at or just below the KK threshold and provides consider- 
able stability. Values of g 2 and the ratio tkk = Qkk/^it are onr y slightly correlated in fitting tjtttt data; there 
is a weak tendency (5%) for them to go up and down together. The data give well defined values for parameters 
of a (980): 

0.9874 ± O.OOlO(stat) ± 0.0030(syst) GeV, (18) 
0.164 ± 0.007 ± 0.010 GeV 2 , (19) 

1.05 ±0.07 ± 0.05. (20) 

These values have changed from earlier publications because of the inclusion of the cusp in ReII(s). The 'edge' 
observed in the t/tttt data due to ao(980) provides a good determination of the coupling to KK. In this respect, 
Crystal Barrel data have an advantage over Kloe data (to be discussed further below). The Kloc data however 
provide an excellent view of the ao(980) line-shape below the KK threshold. A simple program evaluating the 
formulae for ao(980) and the cusp in ReII(s) is available from the author. 

The systematic error on the mass M arises from uncertainty in the mass calibration of the Crystal Barrel 
detector. Systematic errors for g 2 n and tkk arise as described above for the entire fit to r]ir°ir° data. The fitted 
mass of ao(980) is close to the lowest KK threshold, just as the mass of X(3872) is close to the threshold of the 
lowest charge combination in DD*. 

An important detail is that the a (980) could be produced either by its pole term or via the elastic t\tt 
amplitude. Both have been tried, and the fit strongly prefers production via the elastic scattering amplitude, 
i.e. with the Adler zero in the numerator of the production amplitude. 

5 Data for pp -> K^K^ 

The Dalitz plot for these data in liquid hydrogen is shown in Fig. 13. There arc prominent vertical and horizontal 
bands due to if* (890). A detail is that it is necessary to fine-tune the masses and widths of the separate charge 
states for iT*(890). There are also diagonal bands due to a2(1320) and ao(980). The ao(1450) lies in the lower 
left corner of the plot, near the crossing if* (890) bands. As for pp — > t]tt°tt , it is necessary to remove some edge 
bins. There is also a background from pp — ► tt + tt~tt° described in the Crystal Barrel publication. It peaks in 
edge bins close to the left-hand corner of the Dalitz plot. It is necessary to remove 21 bins, leaving 741. 

One would hope to determine the ratio QkkIQt]^ f° r a o(980) from relative contributions in t]tttt and KKir 
data. Unfortunately, when one tries to do this, a serious difficulty appears. It arises from the question of how 
to parametrise the Ktt S-wave amplitude in the KKtt data. This leads to uncertainties in interferences between 
the Ktt S-wave and diagonal bands due to a2(1320), ao(980) and ao(1450). Uncertainties in these interferences 
then lead to uncertanties in the magnitudes of the a 2 and ao signals. To grasp these points, it is necessary to 
review current understanding of the Ktt S-wave, which has advanced a long way since the earlier analysis of the 
Crystal Barrel data [7] . 

In that early analysis, it was assumed that the Ktt S-wave amplitude in production data is identical to that in 
Ktt elastic scattering. Experience with tttt data now makes that appear unlikely [19]. The Crystal Barrel paper 
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Figure 13: Data and fits for pp — > K%K ir^ at rest, (a) Dalitz plot; mass projections for (b) K + tt and K tt + , 
(c) K\tt~ and K\tt + , (d) K^K^ . Points with errors show data and histograms show the fit. 



gives an explicit parametrisation for the K it S-wave. Fig. 14 shows its phase as a function of mass, compared to 
LASS data [32]. In those days, the LASS data were fitted with an effective range expression. Since then, it has 
been recognised that Chiral Symmetry breaking produces an Adler zero in the Ktt S-wave below threshold at 
s = m 2 K — 0.5m 2 . . The dashed curve shows the amplitude fitted to the LASS data including this Adler zero [33]. 
The k pole observed in BES2 [34] and E791 [35] data was fitted simultaneously and is therefore well constrained. 
There is a discrepancy for Ktt masses near 1400 MeV. It is now known that this discrepancy can be removed by 
a full treatment of the cusp at the Krf threshold. For present purposes it is irrelevant since KKtt phase space 
ends at 1381 MeV. 

The discrepancy near the Ktt threshold between the dashed and full curves on Fig. 14 indicates a contribution 
in Crystal Barrel data from the k pole, which peaks near threshold with a half-width of ~ 350 MeV. This is 
confirmed by a fit to the Crystal Barrel data using a 2-component fit to the Ktt S-wave, as in Eq. (8). It 
contains one component from Ktt elastic scattering and a second from the k pole. Uncertainties are compounded 
by the fact that there are contributions from ktt for both 1 = 1 and 1 = initial states, making four Ktt S-wave 
amplitudes in all. These two isospins have opposite relatives signs for coupling to Ktt and are responsible for 
the difference in Ktt distributions between Figs. 13(a) and (b). 

There is in addition the possibility of a third component due to a i^o(1430) amplitude different in magnitude 
and phase in production data and LASS data. That is the case for the E791 data. However, it turns out that 
adding this freedom only improves x 2 by a small amount, ~ 20 and leads to a very ill-defined fit. This extra 
possible freedom is ignored here. 

The earlier Crystal Barrel analysis recognised the need for a low mass Ktt enhancement and parametrised 
it in an ad hoc way which is not consistent with Chiral Symmetry breaking and the Adler zero. In order to 
get an acceptable fit, contributions were introduced from p(1450) and/or p(1700), which can appear in the left- 
hand corner of the Dalitz plot, in the same mass range as the k pole and ao(1450). It is now necessary to try 
to disentangle the complications of this corner of the Dalitz plot. Unfortunately, these complications lead to 
substantial errors in branching ratios. 

Results will be presented first without any contribution from p(1450) and p(1700); then their possible con- 
tributions will be discussed. Either way, there is a large interference between the k pole and the Ktt elastic 
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Figure 14: LASS data for Kit elastic scattering, compared with (i) the Crystal Barrel fit (full curve), (ii) a 
parametrisation including the Adler zero in Kit (dashed curve). 

amplitude. With this freedom, there is great flexibility in what can be fitted to the Ktt S-wave. A free fit 
to a (980) and a 2 (1320) magnitudes gives x 2 = 637.3 for 720 degrees of freedom. If the n pole is omitted, 
X 2 — > 1230; this is clearly unacceptable. The free fit gives a ratio of intensities ao(980)/a2(1320) = 0.0637 com- 
pared with a prediction from fits to rjirir data, Eqs. (18-20), giving a ratio 0.0369 ± 0.0068. This might appear 
to be a large discrepancy. However, if the fit is constrained to agree with the latter prediction, x 2 increases by 
only 2.01 to 639.31, for 2 parameters determining the complex coupling constant A of ao(980). There is clearly 
no disagreement with the prediction; there is just a large flexibility in the amplitudes. 

The same picture emerges for the ratio of intensities for ao(1450)/a2(1320). A free fit to a (1450) gives a 
ratio of intensities ao(980)/a2(1320) = 0.535 compared with its predicted value 0.311. However, again the ratio 
can be fixed to the prediction with a x 2 increase of only 2.01 for 2 less parameters for A1450. 

With the introduction of extra contributions from p(1450) and/or p(1700) the situation is similar. If both 
are introduced, the solution becomes very unstable with excessively large contributions from both resonances 
and destructive interference between them. This destructive interference is a familiar symptom of over-fitting 
the data, so it is necessary to use only one of them. That procedure was adopted in the earlier Crystal Barrel 
analysis. There are recent data from Babar [36] on e+e~ — > jK + K~ir° using Initial State Radiation. They 
observe the same instability between p(1450) and p(1700) contributions, but a dominant p(1450) amplitude. 
Following this lead, the present data are fitted including only p(1450). The story which then emerges runs close 
to that described above. For a fit where ao(1450) and a (980) are constrained to predictions from Eqs. (18-20), 
X 2 drops from 639.31 to 616.50 with the addition of ao(1450). If ao(980) is set free, x 2 improves by 3.09 to 
613.41. This cannot be regarded as a significant improvement. If ao(1450) is set free, the improvement is only 
1.74, again insignificant. 

So the conclusion is that data on pp — > K^K^ir^ are consistent within errors with (i) the parameters of 
a (980) deduced from rjirir data, (ii) the SU(3) prediction for g 2 KK /g 2 ^ of a (1450). However, they do not 
constrain those parameters well. A close inspection of the KK mass projection on Fig. 13(d) shows that data 
favour a slightly narrower ao(980) than is fitted. This could arise from either stronger coupling of ao(980) to 
KK or a steeper form factor, i.e. a larger a parameter and larger radius of interaction. 

A final question is whether there is any significant a (1450) — > KK signal at all. This may be tested 
against the fit where do (1450) and ao(980) are constrained to their predicted values. To err on the pessimistic 
side, this test is made including p(1450) — > KK. Omitting ao(1450) from the fit, x 2 gets worse by 36.8 with 
2 less parameters. This is a change of 5.5 standard deviations, so there does appear to be a signal due to 
do (1450) — + KK. The situation is simply that its magnitude cannot be determined with any precision from 
these data. 
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5.1 Comparison with Kloe data on <fi — > 7(777?") 

There are important data from Kloe on this process giving direct and precise information on the line-shape of 
ao(980). The present status is that data have been published from the first phase of this experiment, with 77 
decaying via both 77 and 7r + 7r~7r° [37]. Preliminary results for parameters of ao(980) have also been presented 
from the second stage of this work [38] . A comparison will be made with these preliminary results. 

The Kloe data define accurately the ao(980) line-shape below the KK threshold. In the process </> — > 700 (980), 
there is a dependence of the cross-section on fc 3 , where fc 7 is the photon momentum. This factor is the usual fc 3 
factor for an El transition. It inflates the lower side of the ao(980) strongly, improving the precision with which 
it can be measured. However, one must beware that the line-shape may be affected by form factors. 

The Kloe group have fitted their data in two ways. The first assumes the KK-loop model of Achasov and 
Ivanchenko [39]. In this model, <j> decays to j(KK), then the KK pair re-scatter via a final state interaction to 
ao (980), which decays to i]tt. With this model, the preliminary parameters for a (980) are 

983 ± 1 MeV (21) 
0.156 ± 0.011 GeV 2 (22) 

1.19 ±0.05. (23) 

These values are quite close to those emerging from Crystal Barrel rirnr data, Eqs. (18-20). 

The alternative fit made by Kloe ignores the constraint of the KK loop model and arrives at a value of g 2 
seriously different: g 2 ^ = 0.096 ± 0.009 GeV 2 . The agreement between the first set of results and Crystal Barrel 
data clearly favours the KK-loop model. The disagreement of the second set illustrates the sensitivity to the 
precise equations used to fit Kloe data. 

The published Kloe data were fitted in earlier work using the KK loop model [40]. That fit has now been 
repeated using the parameters for a (980) reported here. There is a significant contribution to the data (up to 
18%) from (f> — > pir, p — ► 777. From present publications, it is not clear how much of this contribution is eliminated 
by experimental cuts. Excellent fits can be obtained with both the Kloe parameters of Eqs. (21-23) and with 
parameters fitting Crystal Barrel data, Eqs. (18-20) by varying the magnitude and phase of the pn combination. 
The magnitude of the fitted ao (980) signal is proportional to g^g 2 KK and is constrained to reproduce the latest 
branching ratios reported in [38]. The fit is almost indistingushable from that shown in Fig. 3 of [40]. When 
information from the full Dalitz plot becomes available, the pw amplitude can be determined accurately in both 
magnitude and phase. 
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5.2 Pole parameters of ao(980) 

In order to find the pole position of ao(980) it is necessary to parametrise the cusp at the KK threshold in the 
real part of the amplitude. At the cusp, there is a discontinuity in slope, due to the opening of the KK channel. 
If there is no form factor in the KK channel, this cusp can be calculated algebraically with two subtractions at 
the KK threshold. It is given in Rcf. [10], Eqs. (8) and (9): 

ReII(s) ■ Pkk, 1-Pkk . , 2 / 0/1 s 

— =3kk = In— , s>4m KK (24) 

9kk n 1 + Pkk 



tan" 1 !-, s<4m 2 KK , (25) 

S 7T 

where pkk = ^kk I \fs above the KK threshold and v is the modulus of this quantity below threshold; Ukk is 
the momentum in the KK rest frame, and is complex below threshold. 

With the form factor present, an empirical parametrisation is needed, based on Eqs. (24) and (25). An 
excellent fit may be obtained to H(s) from the r/Tr threshold to 1.6 GeV replacing these equations by 

RcII(s) = f{s)j KK + F 2 + F 3 s + F 8 s 2 p (26) 
p 

f[s) = l + F iSr + ksl + F 7 sl (2?) 
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(28) 
(29) 



Here f(s) modulates j with a convergent power series about the KK threshold. [A similar expression is used for 
ao(1450)]. The terms F 2 + F 3 s allow for the double subtraction at the KK threshold. The term F s s 2 is mostly 
concerned with fitting ReII(s) above s = 1.5 GcV 2 . Parameters are tabulated in Table 4. 



Fi 


0.722107 


F 2 


0.052635 


F 3 


-0.230099 


Fi 


-0.638335 


F 5 


1.804376 


F e 


2.165924 


F 7 


1.402098 


F s 


0.172984 



Table 4: Parameters (in units of GeV) fitted to Eqs. (26)-(29). 



The physical region lies at s + ie for all channels, in the limit e — > 0. The tjit, KK and rj'ir cuts may be 

labelled by the signs multiplying i for each channel. The pole closest to the physical region has signs H h 

and lies at M — iT/2 = 989.1 — i40.1 MeV. It is reached from the physical region by going around the end of 
the KK cut and is usually called the second-sheet pole. Tests show that this pole moves little with the KK 
form factor. The strong coupling to KK locks the resonance to the KK threshold and the full width T = 80.2 
MeV is determined by the fitted value of g 2 . If the sign of i for the rj'ir channel is reversed, the pole moves to 
997.4 — i46.3 MeV. The change due to the rj'ir channel is small because the threshold opens at 1093 MeV and 
its coupling is weaker than for KK and rjir. However, its effect is not negligible. 



Sheet Pole M-iY/2 (MeV) 

+ - + (989.1 ± 1.0 ± 3.0) - i(40.1 ± 1.9 ± 2.7) 
+ - - (997.4 ± 1.0 ± 5.0) - i(46.3 ± 1.9 ± 4.2) 
+ + + (920±3±20)-i(93±5±20) 

Table 5: Pole positions on sheets labelled by signs of i in channels rjir, KK and rj'ir. The first errors are statistical 
and the second systematic. 



In Table 5, statistical errors are assigned from corresponding errors for fitted values of M and g 2 . Systematic 
errors arise from (i) neglect of mass differences between the three KK charge combinations, (ii) uncertainty 
about form factors, (iii) isospin mixing with / (980). The first of these is estimated from half the spread of KK 
mass differences. The second is estimated from errors in the exponent of the form factor for the KK channel: 
a = 2.0 ± 0.5 GeV -2 . This estimate is obtained from experience in fitting many sets of Crystal Barrel data for 
wider resonances, where the form factor has a stronger effect. The third error due to isospin violation is unknown 
at present. However, one expects isospin mixing to produce effects small compared with width differences of 
/o(980) and a (980); the second-sheet pole for / (980) has T/2 = 17 ± 4 MeV from current BES II data [16]. 
Systematic effects of the rj'ir channel are hard to estimate without data for that channel. They have been 
estimated as half the difference between the first and second entries of Table 5. 

The position of the pole with + + + signs is further from the physical region than the second-sheet pole. 
It has changed greatly from earlier fits using Breit-Wigner amplitudes without form factors or Adler zeros. In 
that earlier work, this pole, commonly called the third-sheet pole, lay close to 1040 — i83 MeV. The large change 
in its mass arises from sensitivity to form factors. In this respect, Crystal Barrel data fitted here, although 
they determine M and g 2 with modest errors, do not give accurate information on the precise line-shape in rjir. 
This is because of uncertainties in interferences with the irir S-wave and open questions about exactly how to 
parametrise it. Forthcoming Kloe data should improve this situation substantially. The line-shape in those data 
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measures directly the form factor for coupling of rjir to the resonance. The systematic error assigned in Table 

5 assumes the form factor has exponent a = 2.0 ± 0.5 GeV~ 2 , corresponding to a Gaussian source with a root 
mean square radius 0.68 ± 0.08 fm. 

Positions of the poles in the sheets + H — and H corresponding to the rj'n channel will not be given 

because of the absence of data for that channel. 

6 Conclusions 

Experimental conclusions are straightforward. Firstly, dispersive corrections and the s-dependence of amplitudes 
play a major role for both ao(1450) and ao(980). The fit to pp — > ?77r 7r is decisive in settling parameters of 
do (1450). Without these data, the width fitted to a>7r + 7r _ 7r° data alone is unreasonably small ^110 McV. With 
the inclusion of decays of ao(1450) — * ao(980)er, there is an excellent fit to both sets of data, giving the fitted 
parameters of Table 2. The \ 2 01 both fits improve with the inclusion of the s-dependence of the ao(1450) 
amplitude. More importantly, the fit stabilises in a narrow range of parameters and the fitted ao(1450) signal in 
r]Tr°TT° increases by a factor 1.8; this is a clear indication of a better fit to the line-shape. The overall conclusion 
is that ao(1450) decays weakly to r/ir, KK and rj'ir and dominantly to ujp and a (980)cr. It is unfortunate that 
data on pp — > KlK ± n T do not constrain the ratios g'xx/g^ for either ao(1450) or a (980) tightly. The data 
are consistent with the SU(3) predictions. 

For a (980), there is quite good agreement between Crystal Barrel r;7r 7r data and Kloe data. Each have 
their merits. A limitation of present data is the mass resolution of Crystal Barrel near the KK threshold, ±9.5 
MeV. If progress is to be made on isospin mixing between a (980) and /o(980), a mass resolution better than 
0.5 MeV seems desirable, i.e. 10% of the spread of KK masses. 

An incidental result is a better understanding of the dominant J PC = h signal in pp — > (wp)7r. A plausible 
interpretation of the data is presented in terms of a cusp at the cop threshold and a radial recurrence of 7r(1300) 
close to 1730 McV. However, more precise data are needed to clarify this result in cop, pir and [o"7t]l = i channels. 
The Compass experiment could be a good source of such data. Obviously a search for the missing a 's expected 
at higher mass is sorely needed. 

As regards the interpretation of ao(1430), it seems likely to be dominantly an nn state. Black, Fariborz and 
Schechter have pointed out that a (980) may have a radial excitation in the general mass range of ao(1430) [41]. 
Such a radial excitation would almost inevitably mix with the expected tin state. The radial excitation would 
respond to long range meson-meson interactions, and therefore to the oo(980)ct and uip thresholds. These are 
likely to be responsible for pushing the mass of ao(1450) up to that of Ko(1430). 

One of the essential ingredients in fitting all the data is the 7T7t S-wave (and Kir). Empirically it is necessary 
to parametrise it with the 2-component form of Eq. (8). This gives considerably flexibility to the numerator 
N(s) of the amplitude, even though the denominator D(s) is accurately known. Guidance from theory on the 
way the a and k couple in production reactions would be very helpful to experimentalists. 
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